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Abstract 

For a newly discovered iron-based high T c superconducting parent material KFe2Se2, we present an 
effective three-dimensional five-orbital tight-binding model by fitting the band structures. The three t2 g - 
symmetry orbitals of the five Fe 3d orbitals mainly contribute to the electron-like Fermi surface, in agree- 
ment with recent angle-resolved photoemission spectroscopy experiments. To understand the groundstate 
magnetic structure, the two- and three-dimensional dynamical spin susceptibilities within the random phase 
approximation are investigated. It obviously shows a sharp peak at wave vector Q ~ (n, n), indicating the 
magnetic instability of Neel-type antiferromagnetic rather than (n/2, 7r/2)-type antiferromagnetic ordering. 
While along c axis, it exhibits a ferromagnetic coupling between the nearest neighboring FeSe layers. 

PACS numbers: 74.70.Xa,74.20.Pq,71.10.-w 
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I. INTRODUCTION 



Since the discovery of the iron-based superconductor has been reported [111], a series of su- 
perconducting compounds were found, such as 1111 phase (e.g. LaOFeAs yj]), 122 phase (e.g. 
BaFe 2 As 2 11), 111 phase (e.g. NaFeAs [3]), and 11 phase (e.g. FeSe |4] FeTe 150). Recently 
a new superconductor member KFe 2 Se 2 with T c above 30 K (K 8 Fe 2 Se 2 ) [6J] has been reported, 
and attracted considerable interests for its unique insulator and large magnetic moment properties, 
quite different from other iron-based superconducting materials for its parent material KFe 2 Se 2 is 
iso- structural to BaFe 2 As 2 , but chemically close to FeSe. 

The space group of KFe 2 Se 2 at room temperature is I4/mmm (No. 139) |fj]. Its crystal structure 
is composed of edge-sharing FeSe 4 tetrahedra separated by K cations, as shown schematically in 
Fig. 1. 
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FIG. 1. Crystal structure of KFe 2 Se 2 with a tetragonal unit cell (space group I4/mmm), the axes X and Y 
are along the diagonal Fe-Se direction. The crystal structure parameters are adopted from Ref. 



In KFe 2 Se 2 , there are formally 6.5 electrons per Fe, rather than 6 in other iron-based parent 
materials. Thus it is regarded as an electron over-doped 1 1 system. As a consequence, the hole-like 
Fermi surface pockets possessed in other iron-based materials disappear, and only the electron- 
like Fermi surface pockets were observed in recent angle-resolved photoemission spectroscopy 
(ARPES) experiments Jz, 8] and predicted in recent band structure calculations [9, 10]. Thus, the 



Fermi surface nesting, which is widely suggested to be the origin of the striped-antiferromagnetic 
(AFM) (n, 0) spin density wave (SDW) in FeAs-based materials, together with the inter-band 



scattering of the Fermi surface, which is widely believed to be important for the superconducting 
pairing in FeAs-based superconductor, is absent in KFe 2 Se 2 compound. And the strongly three- 
dimensional characteristic of the Fermi surface is also found in KFe 2 Se 2 . 

Early experiment showed that both TlFe 2 Se 2 [11] and TlFeS 2 lll2n are AFM compounds. In a 
recent experiment, the high-temperature magnetic susceptibility decreases with the decrease of the 
temperature, indicating that the Fe spins interact with each other through the AFM coupling [13]. 
However, the magnetic ordering structure is still an debating issue, as addressed in what follows. 
The local density approximation (LDA) calculations suggested that KFe 2 Se 2 is a striped AFM or- 
der, same to the 1111 and 122 phases of the FeAs-based materials [ k]], while some other authors 
|9| thought it to be bi-collinear AFM with (tt/2, tt/2) wave- vector, similar to FeTe [5]. On the other 
hand, for Fe-deficient materials, K .8Fei.6Se 2 etc., the block checkerboard AFM order is observed 
experimentally, in which its origin is attributed to the ordered Fe vacancy Ill4l.ll5n. Moreover, the 



density function theory (DFT) calculations suggested the TlFe 2 Se 2 is a checkerboard AFM order 
Jl^ll . These discrepant results indicate that further investigations on the magnetic ordering struc- 
ture are urgent for understanding the unique properties in AFe 2 Se 2 (A=K, Tl, or Cs) compounds. 

In this paper, we have made our efforts to resolve this debate. Based on our LDA calculated 
results, we propose a three-dimensional five-orbital model for KFe 2 Se 2 , and show its dynamical 
spin susceptibility diverges at wavevector Q=(n, n, 0), suggesting C-type AFM magnetic insta- 
bility, i.e. a Neel-type AFM coupling in the a^-plane and a weak ferromagnetic coupling along 
the c-axis. This paper is organized as follows: the LDA band structure, Fermi surface and a five- 
orbital tight-binding model fitting to the band structures are presented in Sec. II; the dynamical 
spin susceptibility obtained in the random phase approximation (RPA) is shown in Sec. Ill; the 
last section is devoted to the remarks and summary. 



H. TIGHT-BINDING FITTING OF THE BAND STRUCTURES 

The parent material KFe 2 Se 2 has a tetragonal layered structure with Fe atoms forming a square 
lattice in the high temperature phase. The experimental lattice parameters are a=Z?=3.9136A and 



c= 14.0367 A at room temperature [(]]. There are two Fe atoms per unit cell, where each Fe is 
tetrahedrally coordinated by Se. We have performed a band structure calculation of KFe 2 Se 2 us- 
ing the full-potential linearized augmented plane-wave plus local orbitals (FP-LAPW+lo) scheme 



implemented in the WIEN2K package 



The band structure calculations for other iron-based 



materials showed that the electronic structure is sensitive to a small structural distortion. Thus, 
in order to compare with the experiments, we adopted the experimental structural data [6]. The 
density of states (DOS), the Fermi surface and the band structures within the LDA are displayed 
in Fig. 2, Fig. 3 and Fig. 4, respectively. From the DOS in Fig. 2, it is obviously found that 
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FIG. 2. The total and partial density of states of Fe obtained by LDA. 

this compound is a metal in high-temperature phase, and the Fermi surface is mainly contributed 
by three orbitals, XZ, YZ and XY. Meanwhile the d 3Z 2_ R i and d x i-r- bands mainly distribute from 
-2.6 eV to E f , with a little contribution to the Fermi surface, suggesting that the KFe2Se2 is essen- 
tially a five-band system, similar to other iron pnictides. The Fermi surface of KFe 2 Se 2 is plotted 




FIG. 3. Fermi surface obtained within the LDA for parent material KFe2Se2 in folded Brillouin zone. 

in Fig. 3, displaying a three-dimensional characteristic with four electron-like Fermi surface pock- 
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ets at corner, and two electron-like Fermi surface pockets along Y-Z direction in folded Brillouin 
zone. The hole-like Fermi surface pockets at Y are absent, and only the electron-like Fermi surface 



pockets are presented, w 
structure calculations 



hie 



i is in agreement with the recent ARPES experiments [3, [8|] and band 
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The band structures of the undoped KFe 2 Se 2 in high-temperature phase are shown in Fig. 4. 
To implement our further study, we fit the band structures with a five-orbital tight-binding model. 
Note that the orientation of the coordinate system is chosen so that Fe-Fe bonds are directed along 
the x and y axes, in which x and y axes are rotated by 45 degrees from the X-Y axes, where the X 
and Y axes, along the diagonal Fe-Se direction, refer to the original unit cell. The tight-binding 
Hamiltonian for the five-orbital model is described as, 

H ° = Tj ^ C Lr C JP°- niaa ( 1 } 

i,j iacr 
a,fS,<r 

where C] axT creates an electron on site i with orbital a and spin <x, t"j is the hopping integral 
between the i site with a orbital and the j site with f3 orbital, and jj. is the chemical potential 
determined by the electron filling. In the momentum space the Hamiltonian H is expressed as 

H = fcA* + Ta ^))Cl a C ¥(r , (2) 

where r°^(k) is the kinetic energy term, and e a denotes the on-site energy of the a orbital. The 
five-orbital tight-binding fit of the 10 Fe-3d bands obtained by the density functional theory band 
structure is displayed in Fig. 4. 

The model parameters for the five-orbital tight-binding fitting of the KFe 2 Se 2 band structure are 
listed in the following. The on-site energies measured from the Fermi energy for the five orbitals 
are (e lt e 2 , e 3 , e 4 , e 5 )=(-346.7, -346.7, -437.4, -176.4, -927.9), respectively, in units of meV. 
Here orbital indices (1,2,3,4,5) indicate the d xz , d yz , d x 2_f(d X Y), dxy(d X 2-Y 2 ), and J 3z 2_ r 2 components, 



respectively. Similar to BaFe 2 As 2 II 1811 . the hopping integrals along each direction are 



T 11/22 = 2t x ) y cosk x + 2tyj x cosky + 4t^cosk x cosk y 

±2t\\(cos2k x - cos2k y ) + 4t xxy/xyy cos2k x cosk y 

+At xyy i xxy cos2k y cosk x + 4t xxyy cos2k x cos2k y 

+4t[l(cosk x + cosky)cosk z + At\\Xcos2k x - cos2k y )cos(k z ), 



T 33 = 2t 33 (cosk x + cosky) + At 33 cosk x cosk y + 2t 33 x {cos2k x + cos2k y ), 
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FIG. 4. The band structures of the Fe-3d orbitals obtained by the full-potential linearized-augmented- 
plane-wave (LAPW) and its five-orbital tight-binding fitting. The dot lines are the present local-density 
approximation results, the solid lines are the fitting results. The energies are measured from the Fermi 
energy E F =0.409 eV. 

r 44 = 2t x A {cosk x + cosky) + At^cosk x cosk y + 2t A xx {cos2k x + cos2k y ) 
+At 4 xxy (cos2k x cosk y + cos2k y cosk x ) + At xxyy cos2k x cos2k y 
+2? 4 cosk z + 4t xz (cosk x + cosk y )cosk z + &t^ z cosk x cosk y cosk z , 

T 55 = 2t 5 x (cosk x + cosky) + At 5 xy cosk x cosk y + 2t 5 x 5 x (cos2k x + cos2k y ) 
+At 5 xx \,(cos2k x cosk y + cos2k y cosk x ) + At xxyy cos2k x cos2k y 
+2t 5 z 5 cosk z + At 5 xz (cosk x + cosk y )cosk z , 

T n = At l ^sink x sink y + At xx \(sin2k x sink y + sin2k y sink x ) 
+At x 2 xyy sin2k x sin2k y + %t]^ z sink x sinkyCOsk z , 

ji3/23 _ 2it x x sinkyj x + Ait^coskx/ysinky/x 

-Ait x xxy {sin2k y j x cosk x j y - cos2k x / y sink y / x ), 

ji4/24 _ +2it l x 3 sink x / y + Ait l ^cosk y i x sink x i y 

+Ait ] x A xy coskyi x sin2k x iy + Ait xz sink x / y cosk z 
-At 2 x A z sink x iySink z ± Sit l ^ z cosk y / x sink x / y cosk z 
±Sit l x 4 xyz sin2k x / y co sk y / x co sk z - %t 2 x A xy , z sin2k x i y co sk y / x sink z , 
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ji5/25 _ +2it l x 5 sinkyix + 4it l ^cosk x /ySinky/ x 
+ 8 it l ^ z co sk x / y sinkyi x co sk z , 

T 34 = 4t 3 x 4 X} ,{sin2k y sink x - sin2k x sink y ), 

T 35 = 2t 3 x 5 (cosk x - cosky) + 4t 3 xxy {cos2k x cosk y - cos2k y cosk x ). 

T 45 = At 4 ^sink x sink y + At 4 xxyy sin2k x sin2ky 

+2it 45 sink z + 4it xz (cosk x + cosk y )sink z , 

The intra-orbital and inter-orbital hopping parameters up to fifth neighbors of the five-orbital 
model for the fit of the band structure are shown in Table I. 

TABLE I. The intra-orbital t" a and inter-orbital t. hopping parameters up to fifth neighbors of the five 
orbital tight-binding model through fitting the band structures. All the parameters are in units of meV. 
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III. DYNAMICAL SPIN SUSCEPTIBILITY 



In this section we study the dynamical spin susceptibility of the five-orbital tight-binding model 
for KFe 2 Se 2 . The dynamical magnetic susceptibility is calculated for both the non-interaction case 
and the electron-electron interaction case within the RPA. 

1H 



The orbital-dependent dynamical spin susceptibility is given by H19I . 

Xayiq, ico) = r dre ibn < T T S a (q, r) • S r (-q, 0) > (3) 
Jo 

where a and y label the orbital indices, and the spin operator S a (q)=| Y,k,mn C<L(k + q)cr mn C an (k) 
with spin indices m,n, and fi=l/k B T. The physical dynamical spin susceptibility reads ^(q, ico) = 
\ Zffy^ r (q,M, with 

1 _ b«(k)b y ;(k)bl(k + q)b a v *(k + q) 

Xar (q, ico) = -- V • _i_ f (\r i \ FT^ [/( £ v(k + q)) - / (E u (k))]. (4) 

N j-i ioj + E v (k + q) - EJk) 

Here, the matrix elements of the eigenvector b^(k) =< a\/uk > with indices connecting a orbital 
and \x band, are determined by the diagonalization of the tight-binding model Hamiltonian Eq. (2). 

In the presence of the Coulomb interaction, in addition to the kinetic term in Eq. (2), the 
electronic interaction part of the multiorbital Hamiltonian reads, 

H r =U^ n m n ial + U ' ^ n m n m + (U - Jh)^ n iaa n ifi!T 

i,a i i,tr 

a*P a</3 

-j h 2] c) a ^c iai c] n c m + j h 2] c U c li c m c m ( 5 ) 

i i 

where U(U') denotes the intra-(inter-)orbital Coulomb repulsion interaction and J H the Hund's 
rule coupling. Considering the symmetry of the system, we adopt U =U-2J ff . The RPA dynamical 
magnetic susceptibility is described as 

^(q, ico) = X o(q, ico)[I - T Xo (q, ico)V l (6) 

where xo is the bare spin susceptibility defined in Eq. (4), and the nonzero components of the 
matrices Y are given as (T) aa = U, and (T) a *y = J H . Note that the interorbital interaction U 
doesn't contribute to the RPA spin susceptibility [20]. 

To investigate the effect of the k z -dependence of x(q, ico), we study the two-dimensional case 
with fixed k z and q z =0, and the three-dimensional case with integrating over k z and fixed q z , re- 
spectively. Fig. 5 plots the q-dependence of the real part of the bare and RPA dynamical magnetic 
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FIG. 5. Real part of the bare and RPA enhanced susceptibility along the main symmetry directions of the 
Fe/cell Brillouin zone for k z =0 (a), n (b) in the two-dimensional case (k z is fixed and q z =0). Theoretical 
parameters U=0.6 eV, J H =0, or 0.25U (*0), andy8=0.02. 

susceptibilities in the two-dimensional case for k z =0 and n along the main symmetry directions 
of the Fe/cell Brillouin zone, respectively. It is found that the bare spin susceptibility is nearly 
a plateau-like structure, which is almost equivalent both for k z =0 and for k z =n, without obvious 
peak. As a comparison, the RPA dynamical spin susceptibility is considerably enhanced due to the 
presence of the Coulomb interaction, and a small peak appears at Q=(n, n). Furthermore, in the 
presence of both Coulomb interaction and finite Hund's rule coupling, J/z^O, the (n, n) suscepti- 
bility peak is greatly enhanced. Moreover, the RPA spin susceptibility of the k z =n cut is obviously 
larger than that of the k z =0 cut. In addition to a small peak along the (0, 0)-(7r, 0) line, the largest 
peak appears at Q=(tt, tt), indicating the magnetic instability of the iVee/-type AFM order rather 
than that of the bi-collinear AFM order with Q=(n/2, n/2), or of the striped AFM order with Q=(n, 
0). Similar behaviors are also seen in Fig. 6, which displays the contour of the real part of the RPA 
dynamical magnetic susceptibility for fixed value k z =0 (see (a) and (c)), and n (see (b) and (d)) in 
the q x -q y plane in the two-dimensional case. 

Considering the three-dimensional characteristic of Fermi surface in KFe 2 Se 2 , as seen in Fig. 
3, we also present the k z evolution of the magnetic susceptibility. Similar to the two-dimensional 
case, the three-dimensional bare and RPA dynamical susceptibilities are calculated for fixed q z 
values, q z =0 or n. As expected, the peaks of the dynamical magnetic susceptibility are suppressed 
in the three-dimensional case, in comparison with the two-dimensional result. To clearly show 
the magnetic instability in the three-dimensional case, we displays the susceptibilities at U= 0.7 
eV, as shown in Fig. 7. We find that the peak grows up and the divergent tendency of ^(q) 
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FIG. 6. Contour of the real part of the RPA dynamical magnetic susceptibility of k z =0 ((a) and (c)), and 
n ((b) and (d)) in the two-dimensional case (k z is fixed and q z =0). The theoretical parameters: Coulomb 
interaction U=0.6 eV, J#=0 ((a) and (b)), and 0.25U ((c) and (d)), and £=0.02. 

becomes strong with the increase of the electronic correlations. In realistic compound K r Fe 2 Se 2 , 
the Coulomb interaction U is obviously larger than 0.7 eV, implying that the divergency of the 
^f(q) is much stronger than the present situation shown in Fig. 7. Thus the strong divergency in the 
dynamical magnetic susceptibility suggests the magnetic instability or spin ordering in K,Fe2Se 2 . 




FIG. 7. Real part of the bare and RPA dynamical magnetic susceptibility along the main symmetry direc- 
tions of the Fe/cell Brillouin zone for q z =0 (a), n (b) in the three-dimensional case. We adopt U=0.7 eV, 
J#=0, 0.25U (*0), and 0=0.02. 
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The dynamical susceptibility in Fig. 7 shows a small peak in the (0, 0)-(7r, 0) line and a large 
peak around at (n, n), similar to the two-dimensional case. The sharp peak in the spin susceptibility 
at q z =0 indicates a (n, n, 0) three-dimensional magnetic structure with ferromagnetic coupling 
between the nearest-neighbor FeSe layers in KFe 2 Se 2 , i.e. it shows a ferromagnetic coupling along 
the c-axis and an AFM coupling in the ab-plane (C-type AFM). Notice that the interlayer coupling 
is weak due to the long distance between the nearest-neighbor Fe layers along the c-axis. Other 
LDA calculations J^l for both CsFe 2 Se 2 and KFe 2 Se 2 also suggested a ferromagnetic coupling 
along c axis, in agreement with our result. However, the structures of the spin susceptibility of 
q z =0 and n remain qualitatively the same, as shown in Fig. 8. 




FIG. 8. Contour of the real part of the RPA dynamical magnetic susceptibility for q z =0 ((a) and (c)), and 
n ((b) and (d)) in the three-dimensional case. The Coulomb interaction U=0.7 eV, J#=0 (a) and (b), 0.25U 
(c) and (d), and /3=0.02 are adopted. 



On the contrary, the first-principles electronic structure calculations suggested the striped AFM 
in KFe 2 Se 2 [10], or the bi-collinear AFM ordering in AFe 2 Se2_(A=K, Tl, or Cs) [9]. Though 
Zhang and Singh's Neel-type AFM configuration in TlFe 2 Se 2 lllql agrees with our results in the 
a^-plane, the magnetic coupling along the c-axis is different. The origin of the difference in these 
ab initio results is unknown. In our study, we deal with the electronic correlation in the framework 
of the RPA. Considering the fact that the FeSe-based materials with large magnetic moments may 
be intermediate or strong correlated systems, a proper treatment on the electronic correlation in 



11 



these FeSe-based compounds may be important for understanding its magnetic ground state. 

On the other hand, it is found in the recent experiments that the ordered Fe vacancies in 
Ko.8F e i.6Se 2 lead to a block checkerboard AFM order in the square Fe layers with a large magnetic 
moment 3.31 jUg and T N ~ 559 K Jl^j . then an Fe vacancy order-disorder transition occurs at T s 
~ 578 K. In addition, it is also shown that the AFM order is reduced [21] with the increasing 
Fe content. The block checkerboard AFM order observed in the experiments of K0.8Fe1.6Se2 is 
confirmed in recent LDA calculations I22Q. Meanwhile, AFei.5Se 2 (A=K, Tl, Rb, or Cs) displays 
a A-collinear AFM by the LDA calculations li23ll . Obviously, our prediction on the C-type AFM 
magnetic ordering in KFe 2 Se 2 is different from these experimental observations, showing that Fe 
vacancy with the lattice dynamics plays a key role in these Fe-deficient materials. 



IV. SUMMARY 



In summary, we have presented a five-orbital tight-binding model for newly found KFe 2 Se 2 , 
the energy dispersion and Fermi surface are also given. In addition, we also investigated the spin 
fluctuation through the multi-orbital dynamical spin susceptibility within the random phase ap- 
proximation. The results demonstrate a divergent peak appears at the wavevector Q=(n, n, 0), 
indicating a C-type antiferromagnetic ordering in the parent compound KFe 2 Se 2 . Future experi- 
ment for parent materials KFe 2 Se 2 are expected to confirm this kind of magnetic ordering. Further 
theoretical studies on the dynamical magnetic susceptibility for the Fe-deficient compounds are 
also deserved. 
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